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Because of its weak binding the deuteron is thought to be the ideal target from which to obtain 
information about the deep inelastic structure of the neutron. An accurate determination of 
neutron structure functions is vital for testing various parton model sum rules such as the 
Gottfried sum rule, which measures the flavour content of the proton sea, or the Bjorken 
sum rule in spin-dependent scattering. Recent precision experiments on unpolarised ^ and 
polarised deuterium targets have enabled both of these to be tested more rigorously 
than ever before. However, the improved accuracy of the experiments has also meant that 
a definitive determination of these sum rules requires serious consideration of the nuclear 
corrections that arise when extracting neutron structure functions from the deuteron data. 
Considerable attention has already been focussed on nuclear shadowing [0, Q and meson- 
exchange P, ^ effects, which are typically a few percent below x ~ 0.1. On the other hand, 
the most obvious nuclear corrections are those arising from the fact that nucleons in the 
deuteron are bound, and thus off-mass- shell (p^ 7^ M^). Despite this, a fully consistent 
treatment of off-shell effects in the deuteron is still outstanding. 

The discussion of nucleon off-shell effects is usually based on some ansatz relating the deuteron 
structure function to a convolution involving the free nucleon structure function at a shifted 
value of a; or Q^ B H- Instead, we present a systematic treatment of off-shell effects 
starting from the most general Dirac structure of the off-shell nucleon tensor and a covariant 
DNN vertex function. The general formalism was developed in Ref. 0. Here we show 
how to extract the usual convolution from the full calculation — which is exact in impulse 
approximation (c.f. Ref . [PHH ) . The practical importance of the approximations needed to 
arrive at the convolution formula is estimated within a specific model for the structure of the 
nucleon and the DNN vertex. 

As shown by Melnitchouk et al. 0, in the Bjorken limit the most general form of the 
truncated (off-shell) nucleon structure function is xip^) = -^ Xo + / Xi + 4X2, where Xi 
are real functions of p and q. To calculate the complete deuteron structure function one 
must also describe the off-shell nucleon — deuteron interaction, which we write in the form: 
A{P, p) = IAo{P, p) + 7"^ia(P, p) . Taking the trace of A{P^ p) ® xip, 1) ^-nd integrating over 
p gives the spin-averaged quark distribution per nucleon in the deuteron: 

q^{x) = ^ J dy dp^ {Aq xo + -Ai-pxi + A • g X2) (i) 

where x = —q^/2P ■ q is the Bjorken scaling variable, and y = p ■ q/P ■ q is the fraction of 
the light-cone momentum of the deuteron carried by the nucleon. In impulse approximation 
the deuteron structure function involves an on-shell spectator while the off-shell nucleon is 



struck. The DNN vertex of Buck and Gross |TT[ is therefore ideal for this purpose. For a 
deuteron at rest, it relates ^o and Aia (which include the off-shell nucleon propagators) to 
the relativistic deuteron wavef unctions: 

A = 2MD7r^ M {C - 2V + V} 
Ala = 2Md7i^ I b« - ^ ^^ Pa) C (2) 
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where 

C = u'{\p\) + w'{\p\) + ^,2(|p|) + v'A\p\) 

V = v^{\p\) + t;.^(|p|) (3) 
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The functions u and w correspond to the usual non-relativistic 5*- and D-state deuteron wave- 
functions, while Vg and Vt represent the small singlet and triplet P-state wavefunctions which 



are purely relativistic in origin. In our calculation we use the wavefunctions from Ref.|ll 
with pseudoscalar vr exchange, although similar results are obtained with wavefunctions de- 
rived from a model with pseudovector irNN couplings. (The momentum p is calculated as 
in Eq.(|^) below.) 

Of particular interest is the connection between Eq.(|lD and the usual convolution formula 
P, ^ H, in which the nuclear structure function is a one-dimensional convolution of the 
momentum distribution of a bound nucleon in the deuteron {if{y) in Eq.(§) below) with the 
on-shell nucleon structure function {q^ {x/y)). In terms of the truncated nucleon functions 
the latter is defined as: 

g^(x/y) = AM xT + 4M2 xT + 4p ■ g xT (4) 

where XT' are obtained by setting p^ = M^ and p^ = in the fully off-shell Xi- ^^ order 
to derive such a convolution formula one needs to factorise q^{x/y) from the integrand in 
Eq.(|l]). This in turn requires that either the functions Xi are proportional to each other, or 
two of them equal to zero. Alternatively, the x and A parts of Eq.(|l|) could be factored if 
the functions Ai satisfied Ao/M = Ai ■ p/M^ = Ai ■ q /p ■ q. From Eqs.(|) and (||) it is 
evident that this condition is fulfilled in the p"^ = M^ limit by the term proportional to C. 
By further separating x iiito on- and off-shell parts, Eq.(|^) can be recast in the form: 

q^(x) = r ^ ^(y) q^{x/y) + S^^k'^ix) + S^'^^q'^ix) (5) 
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where the first term is the usual convolution result, in which q is given by Eq.(H), and 
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Here the energy Ep is 
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in the deuteron rest frame, and p^ax — V^d ~ yM'^/{l — y) is the maximum kinematic value 
of p^. 

The two correction terms in Eq. (^ can be identified with the off-shell and P-state components 
of the relativistic DNN vertex, 
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and with the off-shell part of the truncated nucleon structure function, 

5Wg^(a:) = ^^jdydp'(A,xl" + A^■px'i" + A.-qxl") 
where x°/^ = Xi- xT- 



(9) 



We should note that the convolution term which we have identified is not unique. An al- 
ternative definition of ^{y) in terms of C — P, with an analogous redefinition of the d'^-^^q^ 
term, can also produce the required proportionality of ^o and Aia needed to satisfy the 
convolution criterion. Nevertheless Eq.(^ is the most natural choice and one that has been 
followed in many phenomenological treatments. In particular, (p{y) includes the infamous 
relativistic flux factor, (Ep/po) y [l^, and therefore satisfies the normalisation condition 

»i 



dy v{y) 
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This condition is necessary in order that (for the valence component) the convolution term 
alone preserves the baryon number of the deuteron. Since both of the correction terms in 
Eq.(^ are proportional to either (p^ — M^) or the P-state wavefunctions, we expect the 
^(-4)g£* Qj^f^ S^^^q^ terms to be small. However, to obtain a quantitative estimate of these 
corrections requires a model for the off-shell nucleon functions Xi- 



Our model for the nucleon structure function is motivated by the work of Refs. |T3|, |Tj, [T^. In 
particular, we suppose that at some relatively low scale {Qo ~ 400 MeV) the nucleon is well 
described in terms of its valence quarks. Then, as explained in Ref. |^, we can calculate the 
Xi in terms of a set of phenomenological vertex functions, (^^^\p, k), describing the process 
A^ — > q{qq), where the (gg)-pair has spin 5* = or 1. On the basis of bag or constituent 
quark models we expect these diquark states to have masses rris — 700 (900) MeV for S = 
(1) |T^. For simplicity we choose a single function for each type of vertex, say /<l>(^) and 



757. 



a^ , and parameterise their momentum dependence in the form 



$(^)(p,fc) = iV(p^) ,,, \,^^„^ . (11) 
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With the above values for ms (note that as the spectator is on- mass-shell, ?n| = (p — k)'^) 
the cut-offs. As, and exponents, ns, are chosen to fit the on-shell data. The function N{p^) 
parameterises the nucleon off-shell dependence of the vertex functions. For the deuteron, 
because of the strong peak in the px distribution at small transverse momentum {pt ~ 25 
MeV), Nij)^) is well approximated by a constant. 

We could equally well choose to parameterise the distributions at larger Qq, however the 
identification of uis with diquark masses would be lost, as would the use of approximate 
SU(4) symmetry to relate the spin-flavour distributions to the valence quark distributions 
(namely d^ai = qlJi ^^^ u^ai = [qoai + 3^'^/) /2, where qlJi and g|,J are the distributions 
with S=0 and S=l diquark spectators, each with its first moment normalised to unity). 
Convergence of the k^ integrations in the functions Xi imposes constraints on the exponents 
ns'- fif) > 1, 111 > 1.5. The correct valence d/u ratio at large x also requires that rii ^ nQ + 1. 
The values rao(i) = 1.3 (2.4) are found to reproduce the large-x data well when evolved to 
Q^ ~ 5 GeV^. The other two parameters, the cut-offs Ao{i), are 750 (450) MeV. The fact 
that Aq < Ai can be understood in a constituent quark picture from the smaller radius of 
the scalar two-quark system. 

For the sea component of the nucleon we take a similar form for the vertex functions, but 
assume that sea quark distributions are associated with intermediate spectator states having 
somewhat larger masses, ~ 27?!™/. (In bag model calculations [^ the sea corresponds roughly 



to four partons in the intermediate state, compared with two for the valence quarks.) To 
suppress the large-fc^ components of the vertex functions for the sea we choose no(i) = 4 
(5), which, with cut-offs Ao(i) = 0.7 (0.7) GeV, give the observed rapid fall-off at large x. 
In addition, we parameterise the observed flavour asymmetry in the proton sea by taking 

dsea - Usea = 2(1 - r)qsea (0 < T < 1), SubjCCt tO the Constraint Jq dx {dsea - Usea) ~ 0.11 



g, where g^ea = {Usea + 4ea)/2. 

Although at the quark model scale, Qo, the low resolution means that valence quarks dom- 
inate, the work of Gliick et al. [^ suggests the phenomenological need of a small amount 
of glue. For the shape of the input gluon distribution we use their valence-like parameter- 
isation: xg{x) ~ x^(l — x)'^. With the above parameters the second moments of {u + d)yai 
and Qsea are ~ 83% and 3%, respectively, which is sufficient to saturate the momentum 
sum rule. The smaller fraction of momentum carried by gluons here compared with the 
analysis of Ref. |T5| reflects the slightly smaller value of Ql (= (0.39 GeV)^) in our fit, and 
the rapid rise of {x)g with Q^. In Fig. 2 we show the resulting proton structure function, 
F2p = X {Au^ai + dyai) /Q + x{6r + A)qsea/Q (with r ^ 0.91), evolved to Q^ = 5 GeV^. 
Clearly the data is very well described over the entire range of x. 

Without introducing any other parameters we can now calculate the deuteron structure 
function, F2D = 5x (wf^^ + rf^„;) /9 + 20xq^^j9. The resuh, evolved to Q^ = 5 GeV^, 
is also shown in Fig. 2. Here gf^^ = (ufea + '^feaj /2, and the flavour distributions q^ are 
defined to be those for a bound proton in the deuteron, with charge symmetry assumed 
for the neutron distributions (c.f. Ref. [Q). The agreement with recent data from SLAG, 
(reanalysed) EMG-NA2 and NMG F2D data ^ is clearly excellent. 

In order to exhibit the effect of binding and Fermi motion, in Fig. 3 we show the ratio of F2D 
to F2N (= F2p + F2n) as a function of the variable xn = 2x. It shows the same characteristic 



features observed in the nuclear EMG effect for heavy nuclei |T8[. The combined nuclear 
effect in deuterium, due to binding, Fermi motion and nucleon off-shellness, is predicted to 
be about 5% at xn — 0.6—0.7, which happens to be similar to that found in Ref. [|19|] . (Reliable 
predictions for the ratio below xn ~ 0.2 would require inclusion of additional mechanisms 
beyond the impulse approximation considered here.) 

Having obtained a good fit for the proton and deuteron structure functions we are now in 
a position to test the numerical importance of the off-shell corrections in Eq.(^. In Fig.4 
we plot the ratios of (the valence components of) the individual terms to the total F2d- 
Both q^{x) and the convolution component in Eq.(^ are normalised to 1 by choosing N{p'^) 
in Eq. ([TT|) to be a constant, about 0.7% (2.6%) larger for the S' = (1) vertex than the 
corresponding normalisation constant for an on-shell nucleon. As a result the first moments 
of (5(-^)g'D (negative) and S^^^q^ (positive) cancel. The ratio of the convolution component to 
the total deviates from unity by about 1-2% for xn < 0.9. (In fact, for xat > 1 the off-shell 



corrections start to dominate.) Although the off-shell effects are numerically small, in any 
precision analysis of the neutron structure function all nuclear effects should be included. 

In summary, we have calculated the complete relativistic deuteron structure function within a 
covariant formalism, including the effects due to binding, Fermi motion and nucleon virtuality. 
The off-shell effects associated with the nucleon structure function as well as with the DNN 
vertex function are found to give corrections to the usual convolution formula. Consequently 
the common deconvolution procedure 0, ^ of extracting the neutron structure function 
from deuterium data will introduce errors in F2„ of ~ 1 — 2% for xn < 0.9. Off-shell effects 
also give rise to corrections to the convolution formula for the polarised deuteron structure 



function qid p[, and the consequences of these for gin are currently under investigation ^1 . 
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Figure Captions. 

1. Deep inelastic scattering from a deuteron in the impulse approximation. The respective 
photon (g), quark (k), nucleon (p) and deuteron (P) momenta are indicated. 

2. Total proton and deuteron F2 structure functions, calculated and evolved to Q"^ = 5 
GeV^. The data are a compilation of the SLAG, re-analysed EMC-NA2 and NMC data [l^ 
for 4 < g2 < 6 GeV^. 

3. Ratio of the deuteron and isoscalar nucleon structure functions in the valence dominated 
region at Q^ = 5 GeV^. 

4. Off-shell nucleon and DNN vertex corrections, 5'^^^q^ (dotted) and S^-^^q^ (dashed), 
resectively, as a ratio to the total q^ (see Eq.(D). The solid curve is (convolution/total) 
—1. Note that even though the first moments of 5^^^q^ and S^-^^q^ are equal in magnitude, 
because these are divided by the x-dependent q^, the areas under the dashed and dotted 
curves are not equal. 
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